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Abstract. We investigate two distinct universality classes for probe particles that 
move stochastically in a one-dimensional driven system. If the random force that drives 
the probe particles is fully generated by the current fluctuations of the driven fluid, 
such as when the probe particles are embedded in a ring, they inherit the dynamical 
exponent of the fluid, which generically is z = 3/2. On the other hand, if the random 
force has a part that is temporally uncorrelated, the resulting motion can be described 
by a dynamical exponent z = 2 as considered in previous work. 

Keywords: driven diffusive systems (theory), stochastic particle dynamics (theory) 
1. Introduction 

The effective interaction between probe particles in a ffuid has been a subject of 
theoretical and experimental studies for some time. If a liquid is far from thermal 
equilibrium then the interactions involving probe particles cannot be described by 
classical thermodynamic approaches such as the potential distribution theorem [1] . Even 
the large scale random motion of a single particle cannot naively be assumed to be 
Brownian motion as temporal correlations between successive collisions of the probe 
particle with the fluid particles cannot safely be ignored. At present too little is known 
in general terms about random motion in non-equilibrium fluids. 

Some progress may be achieved by considering stylized lattice gas models of fluids. 
Such models sometimes are amenable to analytical treatment and may thus give some 
insight into universal features of the stochastic motion and the effective fluid-mediated 
interaction of probe particles. A paradigmatic microscopic lattice gas model is the 
asymmetric simple exclusion process (ASEP) [21 [3], the large-scale dynamics of which 
are governed by the Burgers equation [4]. In particular, in one dimension the motion 
of so-called second class particles is well understood. Such a probe particle moves with 
an average speed given by the collective velocity Vc = dj/dp where j is the steady 



Dynamical scaling for probe particles in a driven fluid 



2 



state current of the lattice gas with density p. The fluctuations of the position around 
its mean are superdiffusive with the dynamical exponent z = 3/2 of the universality 
class of the noisy Kardar-Parisi- Zhang (KPZ) equation [5]. This behaviour is intricately 
linked to the nature of the current fluctuations in the KPZ universality class. Moreover, 
the second class particle is attracted to regions of large density gradient (on molecular 
scale), thus serving as a microscopic marker of shocks [6]. 

When a large system contains two probe particles, the local density gradient in 
the density of the fluid particles induces an effective long-range attractive interaction 
between the probes. The steady-state distribution of the distance x between the two 
second class particles was found to decay as x"^/^ for large a; jTj. In [8] it was suggested 
that this algebraic behaviour generalizes to for other driven diffusive systems with 
a non- universal exponent b. This leads to phase separation for b > 2. For the dynamics 
of two probe particles, initially introduced into a pure system at a finite distance, this 
implies an increase of the mean distance (6 < 2) or the formation of bound state with 
finite mean distance {b > 2). In recent work [9] an attempt was made to determine the 
stochastic motion of two probe particles by a Fokker-Planck approach. This approach 
implies dynamical scaling with z = 2 and yields a differential equation for the scaling 
form of the dynamical probability distribution for the distance of the probe particles. 
In particular, one obtains the growth law for the mean distance {x{t)) oc and the 
growth exponent u as a. function of b. 

Here we argue that the analysis of [9j applies under certain conditions to the 
dynamics of probe particles which has a random part that is uncorrelated in time. 
A simple model for which such dynamics applies is given below. On the other hand, 
for probe particles driven by stochastic forces originating solely from the dynamics of 
the surrounding fluid particles, as in the case of probe particles embedded in a ring, 
the scaling behaviour is determined by the dynamical exponent of the fluid. Namely, 
one has z = 3/2 for generic driven lattice gases rather than the non-generic Edwards- 
Wilkinson exponent z = 2. This picture is theoretically founded below in a more 
detailed consideration of the statistical properties of the underlying current fluctuations 
and supported by dynamical Monte-Carlo simulations. 

The paper is organized as follows. In section [2] we define a generalization of the 
ASEP with probe particles, originally introduced in ^U\. Moreover we introduce a 
variant of the model, with the same transitions, but driven by a different source of 
noise. In section [3] we use dynamical scaling to determine the growth exponent u for 
the driven fluid with z = 3/2. In section H] we analyze the dynamical behaviour of the 
modified model which is governed by the z = 2 dynamical exponent. In section [5] we 
reconsider the results of [9] and present our conclusions. 

2. Models for driven fluids 

Following |10] we define our model on a one dimensional (Id) lattice. A hard core 
repulsion of particles is incorporated by an exclusion interaction such that each lattice 
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site can be occupied by either a positive particle (+), a negative particle (— ), or it may 
be occupied by a probe particle (0). The continuous time stochastic dynamics of the 
model is defined by nearest neighbor jumps where the particles are exchanged with rates 

H > - + with rate 1 + AV 

+ 0^0+ with rate 1 

0-^-0 with rate 1 . (1) 

Here AV is the energy of the initial configuration minus that of the final one, defined 
by the nearest neighbor Ising interaction potential 

V = ~'-J2^i^^+i^ (2) 

i 

where Si = 0, ±1 according to the occupation of site i, and — 1 < e < 1. We consider 
the range —0.8 < e < 1 of the coupling constant. We mention that in the case where 
e < —0.8 (not studied here) the current- density relation of a domain of particles exhibits 
two degenerate maxima. This could lead to a phase separation of the two maximal 
current phases within the domains, which makes the analysis more involved. For a 
discussion of this case we refer to reference 

This model reduces to the KLS model [13] in the case where no probe particles 
are present. Then the steady state has an Ising measure |T2]. This is also expected to be 
the local steady-state measure away from any probe when the density of probes is zero. 
As in most studies of this model we consider equal densities of positive and negative 
particles. In Monte Carlo simulations to be discussed below, one considers sites with 
periodic boundary conditions and random sequential update dynamics. 

In the absence of nearest neighbor interaction, e = 0, the dynamics defined above 
reduces to that of the TASEP with second-class particles which are the probes. The 
steady state of this system is fully known [7]. The statistical weight of all configurations 
with a given segment length n between the two probe particles is proportional to Z„, the 
partition function of the TASEP on an open chain of length n in the maximal current 
phase [15]. This observation has been used to show that the probability to find the 
two probes at a distance n from each other decays for large n as n^^^'^ [7\. In addition, 
this result can be used to estimate the currents of particles which go in and out of the 
segment trapped between the two probe particles. One finds that the outgoing current of 
+ (— ) particles through the right (left) probe takes the form of the steady state current 
of the TASEP. To leading order in l/ra it is given by = ^ ^1 + The opposing 

currents, namely that of — (+) particles incoming through the right (left) probe, take 
the form f"" = i (l + ■ ^^^i* 

we are concerned with, namely N ^ n and 

N ^ 00, this current is well approximated by = |. 

Following one expects similar behaviour for e 7^ 0. Namely, that the current 
of + (— ) particles bypassing the right (left) probe takes the same form as that of the 
current in an open segment of the same length, governed by the same dynamics. This 
current is given by (again to leading order in 1/n) 

jT =Jo(e)(l + 6(e)/n), (3) 
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where 



^°^^) = ^' ^^^^ = 2 2(v + e) ' 



and V = yjz^ + 1- For the relevant values of e one has < 6 < 9/4. Also, similar to 
the e = case, the incoming current is given in the large N limit by 

r = Jo(e). (5) 

Equation has been derived in pLG] , where b is proportional to the coefficient of the 
non-linear term in the KPZ equation. 

The dynamics of two probe particles embedded in an infinitely long ring may 
therefore be rephrased as follows: The left probe particle jumps to the left (thus 
increasing the size of the cluster between the two probes by one lattice unit) after a 
random time determined by the current of the (infinite) environment. The mean of that 
random time is = It jumps to the right (thus decreasing the cluster size by one 

lattice unit) after a random time with mean r°"*(n) = l/j'"(?T,) which depends on the 
size of the cluster. The motion of the right probe particles can be described similarly. 
Therefore the distance n between the two probes is a random process which is increased 
by one unit by the incoming current. This happens after a random time with mean 
r+ = T™/2, independent of the cluster size. The distance is decreased by one unit after 
a cluster size dependent random time with mean T~(n) = T°'^^{n)/2. 

We stress that the incoming and outgoing currents are correlated in time and hence 
do not generate a memoryless Markovian random motion of the probe particles. As 
discussed below the dynamics of the probe particles is controlled in this case by the 
dynamical exponent 2; = 3/2 of the KPZ universality class. 

An interesting modification of the model is given by the following stochastic 
dynamics: The motion of the probe particles that lead to an increase of the cluster 
is taken to be a Poisson process with constant mean waiting time r™. The motion that 
leads to a decrease of the cluster size is unchanged in the sense that it is generated by 
the jump events inside the cluster. Thus while the combined stochastic motion of the 
probe particles is also non-Markovian, the noise has a contribution that is uncorrelated 
in time. Below we shall refer to this dynamics as mixed. 

This mixed dynamics can be described within a similar framework of the model 
defined above. Consider a fluid segment of fluctuating length composed of +/— particles 
denoted by A embedded in an environment of "probe particles" 0. The bulk dynamics 
of the particles within the segment is defined as before by the KLS hopping rules. The 
left boundary dynamics of a cluster A with length n is defined as follows: 

A + A with rate jo{e) 

- A A with rate 1. (6) 
Similarly at the right boundary one takes 

A A - with rate jo{e) 

A + A with rate 1. (7) 
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Note that a fluid segment of lengtli n = 1 can only increase in size. This model is 
considered in more detail in section HI 

3. Probe particles in a KPZ fluid 

In the steady state the large distance distribution of the two probes decays algebraically 
(X with the non- universal exponent b given in (jlj). For b < 2 this leads to an infinite 
mean distance as time tends to infinity. In order to calculate the temporal evolution of 
the asymptotic mean distance we make a scaling ansatz 

P{x,t) = x-H-f^fix/t^/') (8) 

for the large scale probability density of the distance x between the probes at time t. 
In this scaling picture the distance is considered to be continuous and correspondingly 
the notation is changed from n to x. The scaling function f{u) satisfies the boundary 
condition /(O) = finite, and the normalization condition P{x) dx = 1. Here a cutoff 
6 is introduced to prevent divergence at x ^ for 6 > 1. This generalizes the scaling 
ansatz of reference [9] . 

In the range b > 1 the normalization condition yields (3 = while for 6 < 1 one 
finds (3 = {1 — b)/z. From this we calculate the mean distance between the particles, 
{x(t)) = xP{x,t) dx. Inspecting the large t behaviour of this integral one finds 

t'/' b < 1 

{x{t)) ~ <J t(2-b)/^ l<b<2 (9) 

log(t) b = 2 

A + t-(^-2)/- b > 2 

Here A is a non- universal constant. Hence for b < 2 the mean distance increases 
algebraically in time with exponent u = 1/z {b < 1) ot u = {2 — b)/z {1 < b < 2) 
respectively. At 6 = 1,2 there are logarithmic corrections. For b> 2 the two particles 
are bound to each other with a finite average distance while for 1 < 6 < 2 the particles 
are weakly bound. 

Since the motion of the probe particle is driven by the fluctuations of the current 
we argue that the scaling form of the probability distribution for the distance between 
two probes should be determined by the dynamical exponent of the fluid. This leads to 
the identification of z with the KPZ dynamical exponent. Let us briefly elaborate on 
this point. Consider the position of, say, the left probe particle. The rate of change of 
its position is given by the difference between the incoming currents into the segment 
confined by the probes and the outgoing currents from this segment at the left end. 
The incoming current originates from a very large KPZ fluid and hence the variance of 
this quantity grows as t^/^ [I7j. The outgoing current originates from the segment of 
finite length, and is highly correlated with the incoming currents. In the following we 
show that the distribution of the position of the probe particle is determined by the 
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KPZ exponent of the incoming current, z = 3/2. To this end we study numerically the 
fluctuations of the incoming and the outgoing currents at, say, the left end separately. 
We define Xi+it) as the number of steps of the left probe which increase the distance 
between the probes up to time t. These are the moves to the left. Similarly we define 
XL_{t) as the number of length decreasing steps up to time t, i.e. the time integrated 
number of steps to the right. Analogously we define the quantities Xji±{t) for the right 
probe. 

We have calculated the variances and the covariances of the four random variables 
Xl± and Xfi± in simulations. We find that all these quantities increase to leading order 

as 

^4/3 

(see figure [I]). This demonstrates that the dynamical exponent is z = 3/2. 
The scaling ansatz ([8]) predicts for 1 < 6 < 3 that the variance cr^ of the distance 
x(t) = XL^(t) + Xpt+it) — XL_(t) — Xij_(t) between the two probes increases with a 
lower power of t, namely, o"^ oc This requires the cancellation of the leading term 

contribution of a^, which is in agreement with the results of the numerical simulations 
presented in figure [I](a), where the amplitudes of all the variances and covariances appear 
to be the same. 

The mean distance is predicted to grow as t^^"^)/^ for 1 < 6 < 2. This is confirmed 
by Monte Carlo data presented in figure [2] for some values of b in this range. At 6 = 2 
(e = 0.8) the mean distance is predicted to grow logarithmically in time. Our numerical 
simulations for this case show that the mean distance remains very small (of the order 
of a few lattice constants) in the time interval which was studied, and it grows slowly 
in time. 

To complete our analysis we measured the cumulative distribution function 

P{x,t) = Prob(distance between the two probes > x at time t). (10) 
For this quantity, the scaling ansatz ([H]) yields 

P{x, t) = t^^-^^/'-^F{xr'^'), (11) 

where F{y) = f°° ri^^f{ri)dri. Consequently, in the case of 6 > 1, F{y) is expected to 
behave as y"''^^ for small y, whereas 6 < 1 implies F{0) = 1. 

The results for the scaling ansatz for the distribution function are shown in figure [3l 
where the quantity P[x,t)ti^^~^^ is plotted against xt~^/^ for some values of 6 > 1. The 
data collapse shows that in this regime the scaling ansatz is correct with z = 3/2 
and P = 0. On the other hand for < 6 < 1 no good data collapse is found. This 
may be understood by noting that at 6 = the leading quadratic nonlinearity in 
the KPZ equation vanishes and one is left with a quartic nonlinearity. According to 
standard wisdom, power counting shows that such a nonlinear term is irrelevant in the 
renormalization group sense. One expects a dynamical Edwards- Wilkinson exponent 
z = 2, even though the fluid is driven. This results in strong crossover effects for small 
b, making the scaling form (ITT!) with 2; = 3/2 invalid. Measurements of the fluctuations 
of the individual probe particles (not shown in this paper) are found to be consistent 
with z = 2. 
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Figure 1. Measured variances and covariances of tlie quantities Xn^, X^-, 

and Xr^ for (a) e = 0.4 (6 = 1.75) and (b) e = -0.6 {h = 0.75). {AB)^ denotes 
{AB) — {A){B). Note that because of the left-right symmetry some of the possible 
combinations are equal to each other. The simulation parameters are N = 1600, 10000 
realizations and equilibration times (see below) icquii = 14286 (a), tequii = 12500 (b). 



A subtle point of the simulation is the proper choice of the initial condition. We 
simulated the process on a ring of (even) sites with equal {N/2 — 1) number of + and 
— particles and two probes. Initially the two probes were placed on neighbouring sites 
and any configuration of the remaining A^— 2 particles was chosen with equal probability. 
The process ([1]) was then simulated up to some time teqmi with the modification that 
the two probes are constrained to remain nearest neighbours, i.e., these particles hop 
together as if they occupied only a single site. During this "equilibration" the fluid 




Figure 2. Average distance of the two probes as a function of t for some values of 
e. Solid lines indicate the predicted behaviour for large t. The data are shown on a 
log-log scale. The simulation parameters are e — 0.4: N — 1000, icquii = 14286, and 
64000 realizations; e = 0: N ^ 7000, <cquii = 100000, and 64000 realizations; e = -0.3: 
N = 1600, tcquii = 15385, and 24160 reahzations. 



approaches its stationary structure in the vicinity of the probes up to a distance of the 
order tj^^:^^ At this point, defined as t = 0, the restriction for the relative position of the 
probes is released and the distance between the two probes is measured. In our analysis 
we consider t <^ tequii so that the probe particles always move within an equilibrated 
region. 

Clearly, there are also other relevant choices for the initial condition, which describe 
different physical scenarios. For example one may consider an initial state obtained by 
letting the two probes move freely and t = is defined as the time they first meet after 
some equilibration time. This case is not studied here. We believe that such a change 
in the initial condition will not change the scaling exponent, although it may result in 
a different scaling function. In order to test this statement we carried out Monte Carlo 
simulations of the e = case with equal number of + and — particles. The initial 
configuration was picked randomly with uniform probability conditioned on having the 
two probes on neighbouring sites. This initial state is equivalent to the above studied 
one with tequii = 0. Results of the simulations (not presented here) yield the same 
scaling exponent 2; = 3/2 as before, but the scaling function is different from the one in 
figure [3lc. 

For more than two probes our dynamical scaling approach suggests the following 
scaling form for the joint distribution: 

P{Xi, X2,---, Xm) = x\x\ ■ ■ ■ xljt~^^P{Xit~^,X2t~^, ■ ■ ■ , XMt'^) (12) 
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Figure 3. Simulation results for the KPZ fluid with two probes. The cumulative 
distribution function P{x,t) of the distance between the two probes is measured and 
the quantity P{x,t)t3'^^~'^^ is plotted against xt^'^l'^ for some values of the parameter 
e. The simulations were carried out using an efficient multispin coding algorithm. 
The simulation parameters are: (a) N = 1600, tequii = Hill, t = 142,284,569, 1138 
and 163712 realizations; (b) N = 1000, tcquii = 14286, t ^ 91, 183,366, 731 and 64000 
realizations = 64000; (c) N = 7000, toquii = 100000, t = 128, 256, 1024, 2048, 4096, 8192 
and 64000 realizations; (d) N = 1600, tcquii = 15385, t ^ 197, 394, 788, 1575 and 24160 
realizations. 



for M + 1 probes with z = 3/2. Here the formation of bound states for b > 2 leads 
to a condensation transition for a finite density of probes. This phase transition was 
discussed in [TO] . 

4. Probe particles with mixed dynamics 

In this section we consider the model of a cluster with mixed dynamics and study its 
temporal evolution. In particular we are interested in the growth law governing the 
length of the cluster. As is evident from the definition of the model, two processes 
control the dynamics of a cluster of length x: at both ends there is (a) a growth process 
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with rate jo{e), and (b) a length decreasing process with rate 1, see (jMZl)- To analyze 
the dynamics of the cluster length we introduce processes for the individual motion of 
the probes as defined above. The growth process at both ends is a Poisson process and 
hence uncorrelated in time. Thus {Xji+(t)) = {Xj^^it)) — {Xji^it))"^ = jot and similarly 
for Xl+. 

The length decreasing processes, (t) and Xl_ {t) are correlated in time through 
the dynamics of the particles within the cluster. Therefore the process x{t) is also non 
trivially correlated in time. The mean rate of change of the length decreasing process, 
d/dt{Xn-it) + XL-it)) = 2j°"S depends on x 

We start the analysis of x{t) by demonstrating that the dynamical exponent of this 
process is z = 2. As before, we measure the variances and covariances of X^j-{t) and 
XL±{t). We find that for all these quantities the asymptotic growth is linear in t. In 
combining these expressions to calculate o"^ the leading order terms cancel for 6 > 1 
and the growth of cr^ becomes sub-linear, just like in the case of the KPZ fluid. The 
fluctuations of the motion of the probe particles are thus determined by the fluctuations 
of the uncorrelated incoming current (see figure H]). The temporal behaviour of the 
average distance between the probes for various values of e is given in figure [^. The 
measured value of the growth exponent u obtained from equation (Q is consistent with 
dynamical exponent z = 2. 

In characterizing the dynamics of the two probes one may argue that they can 
be considered as two coupled random walkers. The driven fluid in between the probes 
generates an effective interaction, thus reducing the model to a two-particle problem. 
These particles hop away from each-other with a constant rate of (j5]). On the other 
hand the rate at which they hop towards each-other, , depends on the distance x 
between them according to ([3]). This description neglects the fact that the two-particle 
process is not memoryless and leads to a reduced master equation for the interparticle 
distance, which can be transformed into a Fokker-Planck equation in a continuum limit. 
This approach was taken in reference [9] and the resulting scaling forms for the distance 
distribution were calculated. Here we compare this analytical result with simulations 
of the mixed model. In the full range of 6 < 2 we obtain good data collapse with a 
dynamical exponent z = 2, which is consistent with the Fokker-Planck approach, see 
figure El However, For 6 < 1 the scaling function does not seem to agree with the 
Fokker-Planck prediction of reference [9]. We have no explanation for this discrepancy. 

An interesting generalization of the mixed model is obtained by allowing for probe 
particles between the fluctuating boundaries of the system. These extra probes have 
the same dynamics as the particles in the KPZ fluid ([T]) while the evolution at the 
boundaries is unchanged ([MT]). Here the extra probes divide the system into several 
clusters. Numerical simulations (not presented here) show that the evolution of these 
clusters is described by the same dynamical exponent 2; = 2 as for the case of no 
intermediate probes. We measured the quantities X^^, Xji+, Xl_, Xji_ in a simulation 

I Similar data for other values of b are given in reference [9] . In that reference the figure was erroneously 
referred to as corresponding to the KPZ fluid rather than to the mixed dynamics. 
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Figure 4. Measured variances and covariances of tlie quantities X^-, 
and Xr_ in the mixed model for (a) e = 0.4 (6 = 1.75) and (b) e = -0.6 {b = 0.75). 
{AB)c denotes {AB) — {A){B). Note that because of the left-right symmetry some of 
the possible combinations are equal to each other. The averages were taken over 10000 
realizations. 

for intermediate clusters and the result is very similar to those shown in figure |H Thus 
the uncorrelated noise acing on the boundaries dominate the entire system, including 
the inner probes which are not affected directly by it. These results suggest a scaling 
form similar to ( |T2l) with z = 2 for a system with M clusters (M — 1 inner probes). 
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Figure 5. Results of numerical simulations for the average distance {x{t)) between 
the probes for various values of b. Lines are drawn with a slope as expected from 
equation [9l 



5. Conclusions 

In this paper the dynamics of two probe particles embedded in a one dimensional 
driven fluid is studied. The non-equilibrium dynamics of the fluid generates an effective 
interaction between the probes. A simple model is introduced to study this dynamics. 
Within the framework of this model, sufficiently strong attraction between the fluid 
particles leads to a bound state of the probe particles. For less pronounced attraction 
(or repulsion) a weakly bound state is found in which the average distance between 
the probes grows algebraically with time, with a non-universal growth exponent. The 
probes become unbound for sufficiently strong repulsion of the fluid particles, leading to 
dynamics characterized by a universal growth exponent determined by the dynamical 
exponent z of the fluid. 

In previous work by some of the authors of the present paper [9] a similar analysis 
was performed using a Fokker-Planck approach, implying a dynamical exponent z = 2. 
Here we argue that this analysis does not apply to the KPZ fluid which was envisaged 
in In order to clarify the validity of that analysis, two cases which lead to different 
dynamical exponents are considered here. 

In the first case the probe particles are embedded in a KPZ like fluid. Here the 
dynamical exponent which determines the motion of the probes is predicted to be the 
KPZ exponent z = 3/2. This is verified by numerical studies. 
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Figure 6. Simulation results for the mixed model. The cumulative distribution 
function P{x, t) of the distance between the two probes is measured and the quantity 
P{x,t)t^^^~^'> (for b > 1) and P{x,t) (for 6 < 1) is plotted against xt~^^'^ for some 
values of the parameter e. The solid line indicates the prediction of the Fokker-Planck 
equation which is in good agreement with the data collapse for 6 > 1. The simulation 
parameters are: (a) t = 4096,8192,16384,32768,65536,131072,262144,524288; (b) 
t = 46811,93622,187246,374491; (c) t = 43691,87381,174763,349525; (d) t = 
40960, 81920, 163840, 327680; and 50000 reahzations for each case 



In the second case the system is coupled to an external driving force which is 
temporally uncorrelated. This leads to mixed dynamics where both the temporally 
correlated noise generated by the KPZ fluid and the uncorrelated external noise drive 
the motion of the probes. Here we observe numerically that the dynamical exponent is 
z = 2. In this case the analysis of [9] applies in the range b > 1. 

The different dynamical exponents of the two models can be traced back to the 
fluctuations of the incoming and outgoing currents. We observe that in both cases the 
exponent is determined by the incoming current. This result is found to hold when the 
two models are generalized to include any finite number of probe particles. It would be 
interesting to understand how the interplay of the correlated and the uncorrelated noise 
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in the mixed model generates z = 2. Also the crossover behaviour for 6 < 1 that leads 
to poor scaling collapse in the KPZ case is an intriguing open problem. 
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